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LONG-TERM  GOALS 

The  long-term  goal  of  this  research  is  to  develop  an  integrated  model  to  describe  the  generation, 
propagation,  and  transformation  of  large  amplitude  internal  solitary  waves  through  interaction  with 
bottom  topography  and  to  investigate  sea  surface  roughness  associated  with  the  internal  wave  motions. 

OBJECTIVES 

The  objective  of  this  project  is  to  investigate  large  amplitude  internal  wave  dynamics  and  the 
associated  surface  roughness  changes  with  developing  an  integrated  internal  wave  prediction/remote 
sensing  model  based  on  (1)  the  strongly  nonlinear  internal  wave  model  initialized  with  in-situ  data  and 
(2)  the  fully  nonlinear  surface  wave  model  coupled  with  the  surface  current  predictions  of  the  internal 
wave  model. 

APPROACH 

Our  approach  is  to  solve  numerically  the  strongly  nonlinear  internal  wave  model  for  a  two-layer 
system  (Choi  &  Camassa,  1999)  derived  via  an  asymptotic  expansion  method  under  the  sole 
assumption  that  the  waves  are  long  compared  to  the  total  water  depth  with  including  the  leading  order 
non-hydrostatic  corrections.  Once  the  surface  current  field  induced  by  internal  solitary  waves  is 
obtained  from  the  numerical  solution  of  the  strongly  nonlinear  internal  wave  model,  the  evolution  of 
surface  gravity-capillary  waves  is  studied  using  a  surface  wave  prediction  model  that  solves  a  closed 
system  of  nonlinear  evolution  equations  (West  et  al.,  1987;  Choi,  1995)  for  the  free  surface  elevation 
and  the  velocity  potential  evaluated  at  the  free  surface. 

WORK  COMPLETED 

(1)  The  strongly  nonlinear  internal  wave  model  suffering  from  the  Kelvin- Helmholtz  instability 
due  to  a  velocity  jump  induced  by  the  interfacial  displacement  is  regularized  by  changing  the 
dispersive  characteristics  of  short  waves  that  the  original  model  poorly  describes. 

(2)  Comparison  between  the  strongly  and  weakly  nonlinear  internal  long  wave  models  is  made  for 
an  internal  solitary  wave  propagating  over  bottom  topography  and  the  effects  of  strong 
nonlinearity  on  both  the  interfacial  displacement  and  the  induced  horizontal  current  at  the  free 
surface  are  examined. 
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(3)  Large  amplitude  internal  solitary  wave  solutions  are  found  for  a  two-layer  system  where  each 
layer  has  a  constant  buoyancy  frequency  and  are  compared  with  those  for  a  system  of  two 
constant  density  layers. 

(4)  Direct  numerical  simulations  of  surface  gravity  waves  interacting  with  the  surface  current 
induced  by  an  internal  solitary  wave  are  conducted  to  investigate  the  generation  of  short  waves 
crucial  for  microwave  radar  backscatter. 

RESULTS 


1.  Regularization  of  the  strongly  nonlinear  internal  wave  model 


With  the  rigid  lid  assumption  for  the  top  boundary,  the  strongly  nonlinear  model  can  be  written  as 


1ht+(rhui)x=  0’  'll  =  fh  ~  C,  f?2  =  ^2  +  C 
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where  77 fx.t)  is  the  local  thickness  of  the  i-th  layer,  h,  is  the  unperturbd  thickness  of  the  i-th  layer, 
and  Ui(x,t)  is  the  layer-mean  horizontal  velocity.  When  this  model  is  solved  numerically,  it  is  shown 
in  Jo  &  Choi  (2002)  that  this  model  suffers  from  the  Kelvin- Helmholtz  instability  since  u\  is 


different  from  U2.  To  suppress  such  instability,  the  model  is  re-wrtten  in  terms  of  Ui(x,t ) ,  the 
horizontal  velocity  evaluated  at  the  top  and  bottom  boundaries  for  i= 1  and  2,  respectively,  as 
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and  it  is  found  that  the  solitary  wave  solution  of  the  original  model  is  stable  with  the  regulazied 
model,  as  shown  in  figure  1 . 


Figure  1.  Numerical  solution  of  the  regularized  model  for  the  evolution  of  a  single  internal  solitary  wave  of  wave 
amplitude  of  0.8hl  in  a  frame  of  reference  moving  with  the  solitary  wave.  The  initial  condition  is  not  the  exact 
solution  of  the  regularized  model,  but  no  sign  of  the  Kelvin-Helmholtz  instability  is  observed. 
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2.  Comparison  between  strongly  and  weakly  nonlinear  internal  wave  models 

The  comparison  between  the  strongly  and  weakly  nonlinear  internal  wave  models  for  a  single 
internal  solitary  wave  propagating  over  bottom  topogarphy  are  conducted  and  the  results  are  shown 
in  figure  2.  Compared  with  the  strongly  nonlinear  model,  the  weakly  nonlinear  model  overpredicts 
the  displacement  of  the  interface  and  considerable  differences  in  phase  are  observed. 
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(a)  (b) 

Figure  2.  Comparison  of  the  interfacial  displacements  between  the  strongly  (solid  line)  and  weakly 
(dashed  line)  nonlinear  internal  wave  models,  (a)  t=41942;  (b)  t=51920. 

3.  Large  amplitude  internal  solitary  waves  in  a  linearly-stratified  fluid 

By  adopting  an  asymptotic  expansion  technique  similar  to  that  of  Choi  &  Camassa  (1999)  for  a 
system  of  two  constant  density  layers,  approximate  solitary  wave  solutions  of  the  DJL  equation  for 
a  two-layer  system  of  constant  buoyancy  frequencies  (piecewise  linear  density  profiles)  are 
computed  using  a  dynamical  system  approach.  A  result  is  shown  in  figure  3  and  a  detailed 
description  of  such  solutions  for  various  physical  parameters  can  be  found  in  Goullet  &  Choi 
(2008). 


Figure  3.  Streamlines  and  the  horizontal  velocity  at  the  center  for  mode-1  internal  solitary  wave  of  wave  speed 

1/2 

c/(ghl)  =1.2.  The  upper  layer  is  linearly  stratified  while  the  density  of  the  lower  layer  is  constant.  The  density 
across  the  interface  is  continuous  and  the  depth  ratio  of  the  thickness  of  the  lower  layer  to  that  of  the  upper  layer  is  5. 
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4.  Evolution  of  surface  gravity  waves  interacting  with  the  surface  current  induced  by  a 
single  internal  solitary  wave 


without  current:  t=0s 


with  current:  t=0s 


(a) 


Figure  4.  Nonlinear  interaction  between  a  surface  wave  field  chatacterized  by  the  JONSWAP  spectrum  of  peak 
enhancement  factor  y=3.3,  significant  wave  height  #s=0.1  m,  and  peak  wave  period  7p=2  sec,  and  a  surface  current 
generated  by  an  internal  solitary  wave  of  amplitude  a=1.25  m  in  a  two  layer  system  of  h  1=2.5  m  and  /*2=7.5  m.  (a)  Initial 
wavenumber  spectrum;  ( b )  Spectrum  at  t=20s;  (c)  Surface  wave  field  at  t- 20  with  the  surface  current.  The  gereation  of  hi 
wavenumber  components  can  be  seen  clearly  from  the  wave  number  spectrum  at  £=20  with  the  current  as  shown  in  ( b ). 
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IMPACT/APPLICATIONS 


The  project  is  expected  to  provide  a  comprehensive  tool  for  predicting  and  monitoring  the  strongly 
nonlinear  internal  wave  activity  and  surface  signatures. 

RELATED  PROJECTS 

This  investigation  is  closely  related  to  other  projects  funded  by  the  ONR  Non-Linear  Internal  Wave 
Initiative  (NLIWI).  Comparison  of  our  numerical  solutions  with  field  data  for  the  propagation  of 
internal  solitary  waves  and  SAR  images  is  under  progress. 

The  PI  of  this  project  is  working  on  the  University  of  Michigan  MURI  project  entitled  “Optimum 
Vessel  Performance  in  Evolving  Nonlinear  Wave  Fields”  and  investigating  the  energy  transfer  and 
dissipation  of  short  waves  due  to  wave  breaking  whose  understanding  is  important  to  radar  remote 
sensing  of  internal  solitary  waves  is  being  investigated. 
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